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Abstract 

In the context of an Abelian gauge symmetry, spontaneously broken at a 
first-order transition, we discuss the evolution of the phase difference between 
the Higgs fields in colliding bubbles. We show that the effect of dissipation, 
represented by a finite plasma conductivity, is to cause the phases to equlibrate 
on a time-scale, determined by the conductivity, which can be much smaller 
than the bubble radii at the time of collision. Currents induced during the 
phase equilibration generate a magnetic flux, which is determined by the 
initial phase difference. In a three-bubble collision, the fluxes produced by 
each pair of bubbles combine, and a vortex can be formed. We find that, 
under most conditions, the probability of trapping magnetic flux to form a 
vortex is correctly given by the "geodesic rule" . 
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I. INTRODUCTION 



Phase transitions in the early universe, as in condensed-matter systems, may lead to the 
formation of defects, such as cosmic strings To determine the observational implications 
of defects one must follow their evolution from an initial state soon after the phase transition 
up to the time when they begin to have visible effects. There has recently been renewed 
interest in the problem of estimating the initial defect density, which provides the starting 
point for such studies. 

In this paper, we are concerned with the initial number density of strings formed at a 
first-order phase transition proceeding by bubble nucleation. For simplicity, we consider an 
Abelian theory with a scalar Higgs field $, described by the Lagrangian 

L = D^^*D^'^ - ^F^^F"'' - V{^), (1) 

where 1^^$ = 5^$ -|- zey4^$, F^j, = dfj,Ai, — d^Afj,, and the potential y is a function of |$| 
with a local minimum at |$| = and a global minimum at |$| = r]/y/2. 

Nothing happens at the theoretical critical temperature Tc (where the two minima of the 
high-temperature effective potential are degenerate). Once the temperature has fallen well 
below Tc, tunnelling can occur from the false vacuum state at $ = to the true vacuum, 
where $ ~ This is an example of spontaneous symmetry breaking: within each 

bubble there is a random choice of the phase angle 6. 

The bubbles nucleate at random points in space and then expand. The bubble walls 
accelerate, often nearly reaching the speed of light, until they meet and fill the whole of 
space. The nucleation rate per unit space-time volume is determined by the tunnelling 
probability, which may be calculated from the "bounce" solution of the Euclidean field 
theory [Q. It is likely to be a rapidly rising function of time, but for purposes of illustration 
let us assume that it is zero up to some given time and thereafter a constant 7. We also 
suppose that the bubble walls expand at speed v, with 7^ ^ ^ (where H is the Hubble 
parameter). Then the final mean number density of bubbles will be approximately 1/.^^, 
where ^ = (f/7)^/^. (These assumptions are in no way critical to our discussion.) 

The traditional picture of string formation is as follows 0]. Within each bubble the 
phase 9 is essentially constant, but phases in different bubbles are uncorrelated. When two 
bubbles with phases 9i and 62 meet, the sharp discontinuity in the phase in smoothed out to 
become a smooth variation. For energetic reasons, the system will tend to choose the shorter 
of the two paths between 9i and ^2- For example, if 6'i = and < < the phase will 
increase from to 6'2; but if vr < 6*2 < 27r, it will instead decrease from to ^2 — Stt. This is 
often called the "geodesic rule" . It has an obvious extension to higher symmetry groups. 

When three bubbles meet, a string or vortex may be trapped in the region between them. 
Whether this happens or not is determined by the net phase change in going from 9i to 6'2, 
to ^3 and back to 9i. It is easy to see that if the three angles are independent random 
variables, uniformly distributed between and 27r, and if the geodesic rule is correct, then 
the probability of trapping a string is 1/4. Thus the initial string density (length per unit 
volume) is of order 1/4^^. 

This is obviously only a crude estimate. Several of the assumptions we have made are 
questionable. Some of the objections are relatively minor. It is clearly not precisely true 
that the phase within each bubble is constant. There must be random thermal and quantum 
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fluctuations. So the geodesic rule will sometimes break down. When two bubbles meet, there 
may be local random fluctuations that will cause the phase to go the "wrong" way round in 
joining 6i and 62- But the effect of this would merely be a slight increase in the probability 
of trapping a string. Indeed, there could also be a nonzero probability that the net phase 
change would be ±Att, thus trapping a string of winding number ±2. 

As the bubbles grow, closed loops within which strings may be trapped will often be 
formed of more than three bubbles. However, as the bubbles expand further into the trapped 
region, it will generally be split up into several smaller regions. Eventually, the remaining 
gaps will always be flnally closed by the junction of three bubbles. 

There is a particular problem in cases where the velocity of the bubble walls is low. When 
f ^ c we can be assured that no bubble collision can be causally affected by any previous 
collision (because any two points on the expanding light cone are spacelike separated). But 
when u -C c, phase equilibration within a composite bubble might have been completed 
before the two components encounter a third bubble, thus reducing the chances of trapping 
a string. Of course, some strings would still be trapped, because the three coUisions would 
sometimes occur nearly simultaneously, but the string density would certainly be less. One 
very interesting open question (to which we hope to return) is whether in that event the 
relative proportion of small loops to long strings would be changed. At flrst sight one might 
think that if strings are rare then most of them would form small loops If so, that would 
have a dramatic impact on the resulting cosmology. 

Underlying all of these points is a more serious concern: does it make sense to talk of 
phase differences between bubbles in a gauge theory? Indeed by a gauge transformation, the 
phase difference can be set equal to any value we please. It has been argued by Rudaz and 
Srivastava that string formation in gauge theories might be strongly suppressed, although 
a recent analysis by Hindmarsh et al suggests that it is not. 

It is clearly important to understand the process of string formation in gauge theories. 
We begin by studying what happens when two bubbles meet. An important feature of the 
process of phase equilibration (not included in previous analyses) is the role of dissipation. 
We include an ohmic dissipation term in our equations. However, it will be useful flrst to 
consider the problem in the absence of dissipation, especially to clarify some issues relating 
to gauge choice and boundary conditions. 

There are two ways of dealing with the problem of gauge invariance. One is to use only 
gauge-invariant quantities, such as the covariant derivative D^O = 8^9 + eA^. This approach 
provides valuable insights, but it is difficult to apply to the string-formation problem because 
we have to consider non-local variables such as the line integral of D^9. The other is to choose 
a gauge, but that too has its problems. We need a complete gauge choice, not merely a class 
of gauges like the Lorentz gauge within which some gauge freedom still exists. 

There are several gauges we might choose: 

I. The Coulomb gauge. This works well provided we are considering only a limited 
region, outside of which there are no relevant charges. It is problematic in other cases, 
because boundary conditions are needed to solve the Poisson equation for the scalar potential 
= . In practice, this gauge is not as convenient as the others, so we shall not consider 
it further. 

II. The unitary gauge, defined by setting the phase angle 9 to zero. This is very convenient 
within any simply connected region where $ 7^ 0. Unfortunately, it does not work if the 
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region is non-simply connected or contains zeroes of $, i.e., when strings are present. It 
will be useful in discussing the collision of two bubbles, but is less easy to use to study 
three-bubble collisions. 

III. The axial gauge defined hy — 0, together with a choice of gauge on the surface 
z — 0. This is particularly well suited to the problem of two coUiding bubbles, where it 
matches the symmetry of the problem and provides a very neat solution, as we shall see in 
the next section. However, it is not very convenient for more general problems. 

IV. The temporal gauge, given by = 0. This is universally applicable, and very simply 
related to the unitary gauge. To define it completely, we need a gauge choice at some initial 
time, but in our case this is straightforward: we assume that no electromagnetic fields are 
present before bubble nucleation starts, so that we can simply set A = at that time. The 
gauge is then defined for all times. However, this gauge does introduce somewhat artificial 
spatial discontinuities in various functions. 

It will be useful first to examin the collision of two bubbles with no damping of the phase 
oscillations. Although this is a rather artificial problem, the results will help to motivate 
some of the more general discussion later. We begin by defining the gauge-invariant phase 
difference which plays an important role in determining the outcome of the collision. 



II. GAUGE-INVARIANT PHASE DIFFERENCE 

It is convenient to write the Higgs field in polar form: 

$ = -^Xe^'. (2) 

Then the equation of motion for 6 is simply the continuity equation dnj'^ — for the current 

= -eX'D^e = -eX'id^e + eA^). (3) 

The phase gradient D^O is gauge-invariant and can be used to define a gauge-invariant phase 
difference between the two bubbles, 

= dx^DkO, (4) 

where k = 1,2,3 and points A and B are taken in bubble interiors. While the bubbles are 
well-separated, D^9 is appreciably different from zero only in the region where X ?a 0, so 
that and F^^ ki (assuming that F^j, = prior to bubble nucleation). Then 

is independent of the choice of path between A and B. In general of course. A^^ is path- 
dependent, and to make it well-defined, we shall choose the integration path to be along the 
straight line passing through the bubble centers. 

To analyse what happens when the bubbles collide, we shall assume that the radial mode 
is strongly damped, so that X settles rapidly to its equilibrium value rj. The equation for 9 
is then the Klein-Gordon equation 

(a^ + ml)D^e = (5) 
with the gauge boson mass uia — erj. 
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Let us first consider the central region of tlie collision, where the bubble walls can be well 
approximated as parallel approaching planes, so that the problem becomes one-dimensional. 
The initial configuration, just after the bubbles have collided, say at 2; = and t = 0, will 
then have D^O = DyO = and D^O, DtO concentrated near z = within some distance 
5. If the energy of the colliding walls is "instantly" dissipated, then 5 is comparable to the 
Lorentz-contracted wall thickness. (However, if the motion of the walls is relativistic, they 
are likely to overshoot and go through several oscillations before loosing their energy []^ . In 
this case the collision region will be much wider and the phase equilibration process may be 
more complicated). 

It is not difficult to understand that in the course of the following evolution the "wave 
packets" of D^O and Dt9 will spread in the 2;-direction and oscillate at the frequency uo = itla- 
The speed of spreading, u, will be determined by the initial width of the packets: u ~ l/m^5 
for 5 3> and m ~ 1 for 5 ^ ^'a'- (This qualitative picture is confirmed by explicit 
solutions of Eq.(^) with Gaussian "packets"). 

An equation for the phase difference 

M= / dzD.d (6) 

J -ZQ 

can be obtained by integrating Eq.(^ over z, 

(^ + ™i)M^O. (7) 

The boundary terms vanish as long as zq ^ ut. The solution of (^) is A^^ oc cos(m^t). We 
note the "acausal" behavior of A^: it varies on a microscopic timescale ~ m^^ which can 
be much smaller than the radii of the colliding bubbles. 

To extend the analysis beyond the central collision region, we shall now fix the gauge and 
solve the field equations for 9 and A^. Different gauges provide useful, and complementary, 
insights. 



III. COLLISIONS WITHOUT DISSIPATION 

We consider the simplest possible case, where the bubble walls are assumed to move 
essentially with the speed of light. Without loss of generality, we can choose a frame of 
reference in which the two bubbles nucleate simultaneously, say at the points (0,0,0,±i?). 
The bubbles first collide at (i?, 0, 0, 0), when the radii are R. The problem has a high degree 
of symmetry: it is invariant under the 3-dimensional Lorentz group S0(l,2) in the {t,x,y) 
subspace The bubble collision occurs along the surface z = 0, t"^ — x"^ — = R^; for 

any point on that surface there is a frame of reference in which that is the point of first 
contact. 

The obvious gauge to use here is the axial gauge which picks out the one distinguished 
direction, the vector joining the nucleation centres. The symmetry of the problem shows 
that in this gauge we must have 

e{x) = e{T, z), A'^ix) = x°a(r, z), (8) 
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where a = (0, 1, 2) and 

= x'^Xo, = e -x^ - y^. (9) 

The z component of the Maxwell field equation, with the current yields (with a suitable 
choice of arbitrary constant), 

e(T, z) = ArfSa + Tdra). (10) 

Thus there is just one independent unkown function. 

The remaining Maxwell equations, together with the continuity equation for the current, 
show at once that both 6 and A" obey the Klein-Gordon equation with mass erj, whence 
the functions 9 and a satisfy the equations 

d^9 + -d^e - die + e^7]^e = 0, 
r 

d^a + -dra - d^a + e^rfa = 0. (11) 
r 

For the initial conditions we assume that as we go from one bubble into the other, 6 
changes rapidly, say from —Oq to 9o, while 9 = 0. Thus 

^1,=^ = Ooeiz), dj\^=n = 0. (12) 

We also require that A°(x) = up to and on the boundary of the bubble overlap region. 
Using ([T|) and (|10D, this yields 

2 

a\r=R = 0, dra\r=R = —9oe{z). (13) 

rt 

It is straightforward to solve these equations. The solution for 6'a (the subscript indicates 
the gauge choice) is 

6'a = / ^ sin kz (cosuj(t — R) -\ sinu;(r — R)] , (14) 

TXT J k \ ujR J 

where 

= + eW- (15) 
(It is also possible to write this explicitly in terms of the Bessel function of order zero. 



Jo(er/^(r-i?)2-z2).) 

The gauge-invariant phase difference Ad can now be deduced from the asymptotic be- 
havior of ^a(^, X, y, z) at 2 — > ±cxo, 

A9{t) = 9,{t, 0, 0, +oo) - 9^{t, 0, 0, -oo). (16) 

This gives 
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2,R ( 1 ^ 

A6' = — 6*0 cos eriit - i?) H sin eriit - R) 

t \ er]R 



(17) 



Thus phase equihbration does occur in this gauge with a time scale determined, as one might 
expect, by the size R of the coUiding bubbles, with superimposed oscillations with frequency 
given by the gauge-field mass. 
Similarly, Oa is given by 



6*0677^ f dk . 



— sin kz 
k 



T-R 



cosa;(r — R) + i — h 



UJ u 



sinu;(r — R) 



(18) 



It is of course straightforward to transform to other gauges. The unitary-gauge fields are 
given by A';^ = A'^ + {l/e)d^9.^. These fields still exhibit the full symmetry of the problem; 
one finds 



A':, 



x°'b{T,z), b{T,z) 



1 



dla{r,z), 



(19) 



and 



Af, 



--d^e^{T,z). 



(20) 



The unitary-gauge fields again obey the Klein-Gordon equation together with the initial 
conditions 



A 



n\r=R 



0, 



dj3A'^\r=R 



X°'Xf3 2eo 



(21) 



Al\,^R = —^5{z), d^All=R = Q. (22) 

Finally, we can transform to the temporal gauge which, as we shall see, has some rather 
odd features. We find 6^ by integrating the equation 

k = e<. (23) 

We need an initial condition. Assuming as before that outside the bubble overlap region 
6 = 9oe{z), and using (0), we find 

= eoe{z) + e f rbdr = + e f radr. (24) 

JR JR 

It is interesting to note that in spite of our non-covariant choice of gauge, 6t is still a function 
only of r and z. It satisfes exactly the same initial conditions (p!2D as 6'a. However it does 
not satisfy the Klein-Gordon equation. If we integrate the /x = Maxwell equation (the 
Gauss law) using the given initial conditions, we find 

dkA': = er]'[9,-9oeiz)]. (25) 

Thus 9t actually obeys the equation 
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(26) 



The symmetry of the problem also shows up, even more remarkably perhaps, in the fact 
that not only does ^4^ = but in fact — {). The only non-zero component of the gauge 
potential is 



= r Td.adT. (27) 

J R 



Both Al and its derivatives vanish on the initial surface t — R. It is non-zero only because 
it too does not satisfy the Klein-Gordon equation; rather 

d'^Al + eW - ^5(^) j = 0. (28) 

One effect of the extra terms in the field equations is that at late times, 9t, unlike 9s, 
does not become uniform; it tends to the static solution 

lim^t = M-2)(l-e~"'''''). (29) 

t— »oo 

Similarly, 

lim ^Ty^oe"""!^!. (30) 

In a sense phase equilibration never occurs in this gauge. Of course, -D^^t does vanish. 

Clearly we have to be careful if we use the temporal gauge to allow for these extra 
terms in the equations of motion. Note that the time derivatives of 9^ and Al do obey the 
Klein-Gordon equation. 

It is also interesting to examine the form of the gauge fields. We easily find that the only 
non- vanishing components of F^^ are 

F'^' ^x'^d,a{T,z). (31) 

In other words, we have a longitudinal electric field, — —td^a and azimuthal magnetic 
field B"^ = pdztt (in cylindrical polars). 



IV. EQUATIONS WITH DISSIPATION 

We now add an extra term to Maxwell's equations: 

dxF^^ = r + 3^, (32) 

where is the conduction current. This represents the effect of dissipation due to couphng 
with other fields. We shall assume that it can be written in the standard ohmic form, i.e., 
that the spatial part of is given by 

jc = <7E, (33) 
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for some suitably defined conductivity cr. The corresponding cliarge density pc is fixed by 
the requirement that the continuity equation, 



Pc + V-je = 0, (34) 

hold (together with the intial condition pc = at early times). In phase transitions with 
weak supercooling, bubbles expand on a background of a dense relativistic plasma which 
is highly conductive. Even in the case of strong supercooling, the energy released in the 
collision region creates a plasma of its own, and again the conductivity is expected to be 
high. We shall see that phase equilibration between colliding bubbles is strongly affected by 
plasma dissipation. 

Setting A'^ = (0, A), we may write the equations (|32|) as 



-V'^(j)-V ■ A = -eX\e + e(f)) + p, (35) 

and 

A - V^A + VV ■ A + V(/> = eX\ve - eA) - a{A + V0). (36) 

The first of these only serves to define pc and is therefore not very useful. 

We should really also include a damping term in the equation for X. However, we shall 
not consider that explicitly. Instead, we assume that the oscillations in the radial direction 
in field space are efficiently damped so that when the bubbles collide X rapidly settles down 
to its equilibrium value, X = rj. 

Let us then specialize to the region inside the bubbles where X = t]. In this case, it 
seems best to consider first the unitary gauge. (For the moment we drop the subscript u.) 
The equations then become 

+ V ■ A = (37) 

and 

A - V^A + (tA + eVA = -aV(f). (38) 

If the fields vanish outside a finite volume, we can unambiguously separate the longitudinal 
and transverse parts of A, writing 

1 



A = A^ + V-^V ■ A. (39) 

The transverse part satisfies the damped Klein-Gordon equation with mass er]: 

AT - v^AT + aA^ + e^A^ = 0. (40) 
If we look for modes proportional to 

gPt+^k-x^ (41) 

we find for p the dispersion equation 
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2\ 1/2 



P=~±^[k^ + eY-^j , (42) 

representing a damped oscillation. 

For the longitudinal modes, we find the third-order equation 

- + a0 + e^ri^ - aV^ = 0, (43) 

which yields 

+ + e'^7]'^)p + a{p'^ + k^) = 0. (44) 
If a is small, there is one real root, 

P = Pi~ ^, (45) 



where 



and a complex conjugate pair. 



UJ 



a;2 



2 - k^ + eY, (46) 



P = P± ~ - (47) 

For small k, the latter two modes behave exactly like the transverse modes, but the first 
one is a much more slowly decaying mode with no oscillation. 

On the other hand, if a ^ erj and k <^ e^rf /a, then all three roots are real: 

ak"^ e^rf ak"^ e^rf 
Pi ~ — P2 ~ + Pz^-(J + • (48) 

It should be noted that the slowly-decaying mode p\ becomes pure gauge in the limit — > 0. 

In the next section we shall see that this mode does not affect the evolution of the gauge- 
invariant phase difference A^. 

Now let us briefly consider the equations in the temporal gauge. With = 0, the 
equations become 

- V\ - eV ■ At = (49) 

and 

At - V'At + VV • At + a At + e^rfA, = -erfVO,. (50) 

Of course the equation for the transverse part of At is exactly the same as before, eq. (^OD, 
but this time if we eliminate V ■ At using (|49|) we find for 9^ a fourth-order equation: in fact 
9t satisfies the same equation (|^) as does 0u- This is actually obvious because the gauge 
transformation between the two gauges shows that 6^ = e0u- 
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In the temporal gauge we always have the freedom of making an arbitrary time- 
indenpendent gauge transformation. Thus the equations have static solutions in which 9 
is any arbitrary function; of course A must be chosen to make — eA = 0. Note that 
since only 6 appears in the equation, the extra terms discussed in the preceding section do 
not appear; they would reappear, however, if we tried to write an equation for 6 itself, in- 
corporating the chosen initial conditions. In a non-simply connected region, the 6t equation 
allows more general types of solution, in which it is not single- valued. For example, we could 
mimic a closed loop by imposing periodic boundary conditions, and look for solutions with 
9{z + L) = 9{z) + 2n7c. One such is ^ = 2n7rz/L, with = 27111 /eL. This is similar to the 
type of solution one would expect to see eventually around a string, once the oscillations 
have died down. 



V. COLLISIONS BETWEEN PLANAR WALLS 

As in Section II, we shall first consider the central region of a two-bubble collision, where 
the problem reduces to one-dimensional form. 

In one dimension of course the only modes are longitudinal, so we should use ( ^3D or 
its equivalent in the temporal gauge. We have assumed that X settles down rapidly to its 
equilibrium value. It is reasonable then to take as our initial condition, just after the bubbles 
have collided, say atz = Oatt = 0,a configuration where X = rj and (in the temporal or 



axial gauge) 9 has a discontinuous jump from — 6'o to 6'o, as in (O). We also assume that at 



that time there are no electromagnetic fields present, so we may take 

At(0,z) = O, At(0,2) = 0. (51) 

We can now transform back to the unitary gauge, where we find as before the initial 
conditions 

no 

0(O,;^) = O, 0(O,z) = -^5'(^), (52) 



29 

A'U),z) = -5{z), i^ = 0. 

e 



The solution is immediate: 



(53) 



(54) 



where pi and p± are given by (^5|) and (147|) and the amplitudes 921, ip± satisfy 



Piipi + + P-^- = 2ik9o/ e. 



(55) 



The solution is 
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2ikeQ P++P- ._„x 
V'l = 7 V' (56) 

e {p+ - Pi){p- - Pi) 



with two similar equations obtained by permutation of (1, +, - 
The gauge-invariant phase difference (0) can be written as 



oo roo 



oo J — oo 



A^^ / D,ddz = -e A^dz. (57) 



and can be evaluated using Eg. (|5^) . Alternatively, we can obtain an equation for A6'(t) by 
integrating (|36D, 

The boundary conditions corresponding to (^3]) are 

Ae(0) = 2^0, A^(0) = 0, (59) 



and the solution is 



where 



M{t) = A+eP+' + A_eP-\ (60) 



and 



P± = -2±MeV-^l (61) 



A, = T J^. (62) 

P+-P- 



Note that the values of p± are the same as in Eg. (^2]) for the transverse modes at k = 0. 
In the small-o" case, we find 



A9(t) = 2^06""*/^ (coser^t + — siner/t ) . 



(63) 



It is interesting to compare this with the corresponding expression, (|T^), in the undamped 
case. The oscillating factor is essentially the same in both cases. The difference is that 
here the dissipation causes exponential damping, while in the previous case the spherical 
geometry led to a linear damping. For large a there is no oscillation, and Egs.(|60| ) - (|62|) give 

Ae{t) = 2eoexp{-e^r]H/a). (64) 

The characteristic time scale for the current damping, t^, is given by 

td~l/a {a<er]), td-^cr/er] {a>e7]). (65) 
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It is also interesting to look at the solution directly in the temporal gauge. There, we 
have to solve a fourth-order equation, so in addition to (|l^) we need initial conditions for 6 
and 9. These are given by (|5ll) together with (^91): 



6'(0, z) = 29o6'{z), 9{0, z) = 0. (66) 

Thus we find 

9,{t, z)= I ^e'"' (^0 + Ae^'' + ^+e^+* + ^^.e^"*) . (67) 
It is easy to verify, however, that in fact vanishes; for small a, we have 

It is interesting to note that in the limit /c — > 0, — > while i^i ~ — 219^1 k. This shows 
that as before 

lim ^t = ±^o, (69) 

although 9t{t,z) approaches its limit more slowly at late times. 

We can also find a solution corresponding to an idealized one-dimensional version of a 
three-bubble collision, by imposing periodic boundary conditions. We may represent the 
three bubbles at the moment when they meet by the segments [0, Li], [Li, L2] and [L2, L] of 
the z axis. If we choose the initial phases to be 0, 9i,92, where 9i, 92 — 9i and 27c — ^2 all he 
between and vr, then it is natural to impose the boundary condition 9t{L) = 9t{0) + 27r. 
In the unitary gauge, the initial conditions for (p, as in (|52D, become 



0(0, .) = - L,) + ^^^^5'(z - L2) + ^-^^^^6'iz). (70) 

e e e 



The solution is similar to (|5^) but with the Fourier integral replaced by a sum over the 
modes fc„ = 2n7c/L, and the right hand side of ( p5D replaced by a sum of three terms with 
appropriate phases. An interesting special case is the one where Li = L/3, L2 = 2L/3, 
while 9i = 2it/3 and ^2 = 47r/3. Then only modes with n = (mod 3) will contribute. 
The non- vanishing coefficients v^o,3n and v?±,3n are of essentially the same form as in (p6D. 
Note that there is no n = term. The boundary condition on 9t must be accommodated 
by adding to the Fourier series an extra constant term, for example the linear term 2ttz/L. 



VI. MAGNETIC FLUX 

Let us now move away from the immediate neighbourhood of the point of first contact 
between the bubbles. In particular, we want to examine the magnetic fields present in the 
region where the collision occurs. 

For simplicity, we again suppose that the two colliding bubbles have equal radii, R, at 
the moment of collision. The main effect of the curvature is to delay the collision at points 
off the line of centres. If the collision occurs at the origin at time t = 0, then at a transverse 



13 



distance x, it will occur (to first order in the curvature) at time t = x'^/2Rv, where v is the 
velocity of the bubble walls. As a first approximation, we may expect that the fields near 
X = y = have the same form as in the planar case, while more generally, for example, 

e{t, X, y,z)^e(^t- ^^ll^ , 0, 0, . (71) 

This of course has the effect of introducing transverse gradients into the fields. Consequently, 
and jy are non-zero and in turn generate transverse components of A. The physical result 

is the appearance (as in the undamped case) of an azimuthal magnetic field. The magnetic 

flux actually forms a loop around the region where the collision is occurring. 

We are particularly interested in the total magnetic flux. It is possible to calculate that 

without going into the details, by looking at the line integral of the vector potential. 
Consider a rectangular loop ABCD, where (see Figure 1) 

A,B={0,0,tZ), C,D = {X,0,±Z). (72) 

The values of X and Z are chosen so that during the time interval of interest, all four points 
are within the bubbles, but the path CD passes outside them. We have assumed that there 
are no fields outside the bubbles. For simplicity, let us suppose that $ remains slightly non- 
zero, so that the phase 6 is well defined and interpolates smoothly between its values inside 
the bubbles. (It would of course make no difference if we assumed some other distribution, 
so long as $ remains very small, but the argument is simpler in this case.) It is easy to see 
that initially, when the bubbles have just met, 

Dkedx^ = 0. (73) 

ABCD 

In fact, using the definition (0), one sees that A9ab = AOdc = 26*0, while /S.Obc = AO da = 
0. On the other hand, after a time t ^ td, with td given by (|65l), dissipation has reduced 
A9ab to zero, while leaving the other three legs unaltered. Thus 

DkO dx'' = -200. (74) 

ABCD 

By hypothesis the loop contains no zeros of $, so is single- valued. Thus it follows at once 
that 

dz dxBy=i Akdx^ = —. (75) 

-Z Jo JABCD e 

This is exactly what we should expect. Note that when three bubbles meet the three 
magnetic flux tubes combine. If there is a total net phase change of 2tt around the loop, the 
total flux trapped is exactly one flux quantum, 2ii/e. 

It is interesting in passing to see how the result ( ffil) emerges in different gauges. In the 
unitary gauge, of course, each AO is given by a line integral of A; the essential effect is that 
along AB, A^ is non-zero initially, but zero finally. In the temporal gauge, the situation is 
reversed. There is no change in the values of 6 at the four corners. Initially, A^ vanishes along 
AB, but after the dissipation its line integral exactly cancels Ob — Oa- The most complicated 
situation occurs in the axial gauge. Here we find that finally 9a = Ob = 0, which ensures 
that AOab vanishes. However, to keep the line integrals along the transverse sections BC 
and DA zero, A^ must acquire a non-zero value to compensate for the transverse gradient 
of 9; it is easy to check that it does. 
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VII. FLUX SPREADING 



At any particular time t after the initial collision, the magnetic flux forms a ring. We 
can readily estimate its size. The outer edge occurs where the bubbles are just meeting, 
namely at a transverse distance 

X = V2Rvt. (76) 

The time taken for the currents generated to be dissipated at any given point is td given by 
Eg. (|65D . The inner edge of the ring occurs where dissipation is essentially over, namely at 

X = pRv{t-U). (77) 

Note that the ring becomes narrower as it expands. This is because it is slowing down. 
The speed of expansion is \J Rv/2t. If f = c, this is always superluminal (the small-angle 
approximation breaks down of course when t becomes of order R). In that regime, one 
should not think of a tube of magnetic flux moving outwards, but rather of magnetic flux 
being created and destroyed at successively larger radii. On the other hand, if f -C c, the 
speed may eventually be less than c. At that point, we can no longer disregard the region 
outside the bubbles. 

While the magnetic field is expelled from bubble interiors (as in the Meissner effect), it 
has a tendency to spread at the speed of light in the false vacuum region outside the bubbles. 
If the conductivity due to particles present in that region can be neglected, then the outer 
radius of the magnetic flux ring is x ~ t, while the inner radius is given by Eq.(^). At 
t > R most of the flux will escape to distances greater than the bubble radii. If a third 
bubble arrives after that, the probability of vortex formation in this three-bubble collision 
will be strongly suppressed. 

Realistically, however, the case of low bubble velocities corresponds to weak supercooling, 
when bubbles expand on a background of relatively dense plasma, and the conductivity in 
the exterior region is far from being negligible. The magnetic field dynamics in a conducting 
medium is described by the equation 

B - V^B + (xB = 0. (78) 

It is easily seen that the characteristic time of magnetic field variation on a length scale L 
is t ~ (tL^, and thus the speed of magnetic flux spreading is 

Vb ~ L/t ~ iaL)-\ (79) 

In the limit of very high conductivity, vb ^ (the magnetic flux is "frozen" into the plasma). 
If the typical bubble radius at collision is ^ and the velocity is v, then the condition for the 
flux spreading to be negligible is 

> 1. (80) 

In order to assess the validity of this condition, we shall need estimates for C,, cr and v. 
The radius ^ can be estimated as ,^ ~ (17/7)^/"^, where 7 is the bubble nucleation rate. To 
estimate the conductivity cr, we represent the drift velocity of charged particles as 
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■^d ~ p/e ~ eEr/T, 



(81) 



where p is the average momentum gained by a particle between colhsions, E is the electric 
field, r is the mean free time, T is the temperature, and e ~ T is the average particle energy. 
The mean free time is r ~ (n5")~^, where a ~ e^/T^ is the characteristic scattering cross- 
section and n is the density of charged particles. Now, the current density is j ~ nevd ~ 
{T/e^)E, and thus 

a ~ T/e^ (82) 
and the condition (pOD can be rewritten as 

V > e^^T. (83) 

Since the average bubble separation ^ is expected to be greater than the thermal particle 
wavelength T~^, the right-hand side of (^) does not exceed 10~^. The temperature at 
the phase transition is typically T ~ 77, and Eq.(|8^ corresponds to the high-conductivity 
regime, a ^ er]. 

The bubble velocity v in the overdamped regime can be roughly estimated from / ~ Ap„, 
where / is the force of friction due to particle scattering off the bubble wall, and Ap^ is the 
vacuum energy difference between the two sides of the wall. Both / and Ap^ are model- 
dependent, but in all cases / < T'^v, which gives a lower bound on v, 

V > ApjT\ (84) 

Typically, one finds that bubble velocities are rather high. For example, a detailed analysis 
for the case of electroweak phase transition gives v ~ 0.1 — 1. This suggests that the condition 
(^) is typically satisfied, so that the magnetic fiux does not escape. Exceptions may occur 
in cases of nearly degenerate vacua, when Ap„ is very small. 



VIII. CONCLUSIONS 

Our main conclusions can be formulated in terms of the gauge-invariant phase difference 
defined in Eq.(^). We found that A6 remains equal to its initial value while the bubbles 
are well-separated and undergoes damped oscillations after the collision. The oscillation 
period is determined by the gauge boson mass, and the damping time by the plasma con- 
ductivity. Both timescales are typically much smaller than the bubble radii at the time 
of collision. This is to be contrasted to the case of a global symmetry breaking, when the 
phases in the bubble interiors can change only after the arrival of the Goldstone boson waves 
propagating at the speed of light from the collision region. 

The currents associated with the oscillating phase gradients generate a magnetic field 
which is concentrated in a ring-shaped region near the circle on which the bubble walls 
intersect. The total magnetic fiux is determined by the initial value of A6, 

$B = AOo/e. (85) 

Right after the collision, the radius of the intersection circle grows faster than the speed of 
light, but if the bubble expansion is substantially slowed down by interaction with ambient 
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plasma, the speed may eventually become nomelativistic. The magnetic field then begins 
to spread into the exterior region and can in principle escape to distances large compared to 
the bubble radii. However, magnetic field has a tendency to be "frozen" into a conducting 
plasma, and the speed of spreading is typically smaller than the bubble expansion speed 
(with the exception of cases of nearly degenerate vacua, when the bubble velocities are 
extremely small). 

In a three-bubble collision, the three magnetic fiux tubes combine when the bubbles 
meet. The resulting magnetic fiux is $b = ^^'tot/e, where A6'tot = 2?™ is the sum of the 
initial phase differences between the bubbles and n is an integer. For n 7^ 0, a vortex is 
formed carrying n flux quanta, $b = 27™/e. Thus, although the physics of bubble collisions 
turned out to be interesting and complicated, our conclusions are essentially the same as 
those resulting from the naive analysis in Ref . . 

In conclusion, we would like to mention some open questions. As a matter of principle, 
it would be interesting to understand what happens when the bubble expansion speed is 
so low that most of the magnetic flux does escape. One expects that the vortex density in 
this case is suppressed, and we are currently investigating what effect this may have on the 
properties of the resulting string distriution. 

Throughout the paper we assumed that the radial mode of the Higgs fleld is strongly 
damped and settles to its equilibrium value on a timescale short compared to the phase 
equilibration process. In general, however, the energy dissipation in the radial mode is 
itself a complicated process and can take considerable time, especially in the case when the 
colliding walls are relativistic. The phase equilibration in this case can be studied using 
numerical simulations. 

Another interesting and important question is what effect thermal fluctuations of the 
magnetic fleld and of winding number have on the vortex formation. A naive estimate of 
the magnetic flux fluctuation on scale is 5$^ ~ (^T)^/^, where T is the temperature. This 
could in principle be much greater than a single flux quantum. If all this flux is squeezed 
into a vortex by advancing bubble walls, then we can expect that a sizable proportion of 
vortices will have large winding numbers. A particularly interesting situation may arise in 
models where the gauge bosons are heavier than the Higgs particles, so that n > 1 strings 
are stable. In this case one expects the formation of an interconnected network of strings 
with different winding numbers. The cosmological evolution of such a network may be very 
different from the standard picture. 

Thermal fluctuations may also play an important role in defect formation during flrst- 
order phase transitions with a global symmetry breaking. In this case, the phase 6' is a 
massless Goldstone fleld, and its fluctuation on scale ^ is (56* ~ (T^)^/^. Once again, this can 
be large, resulting in the formation of large vortex clusters. We note, however, that multiply 
charged vortices are unstable in the case of a global symmetry breaking, and we expect the 
clusters to disperse. 
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FIGURES 
FIG. 1. 

Two colliding bubbles, showing the closed path ABCD. 
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